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Abstract 

A number of transition-metal (TM) doped group-IV semiconductors, HxMi-x (R=Cr, Mn and 
Fe; M=Si, Ge), have been studied by the first principles calculations. The obtained results show 
that antiferromagnetic (AFM) order is energetically more favored than ferromagnetic (FM) order in 
Cr-doped Ge and Si with x=0. 03125 and 0.0625. In 6.25% Fe-doped Ge, FM interaction dominates 
in all range of the R-R distances while for Fe-doped Ge at 3.125% and Fe-doped Si at both 
concentrations of 3.125% and 6.25%, only in a short R-R range can the FM states exist. In the 
Mn-doped case, the RKKY-like mechanism seems to be suitable for the Ge host matrix, while for 
the Mn-doped Si, the short-range AFM interaction competes with the long-range FM interaction. 
The different origin of the magnetic orders in these diluted magnetic semiconductors (DMSs) makes 
the microscopic mechanism of the ferromagnetism in the DMSs more complex and attractive. 

PACS numbers: 75.50.Pp, 71.55.Cn 
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I. INTRODUCTION 



Diluted magnetic semiconductors (DMSs) have stimulated a great deal of interests be- 
cause of its potential applications in the spintronics, in which the electron spin becomes 
another degree of freedom in addition to the usual charge one. Since the discovery of FM 
order in Mn doped III-V semiconductors such as InAsi and GaAs^^ the transition metal 
magnetic impurities have been used as spin injector, which are doped into the semiconductor 
hosts to cause the ferromagnetism. But usually, the Curie temperatures of these materials 
are far below the room temperature, which, for example, is 35 K for Ini_a,Mna.As and 110 
K for Gai_^Mn^.As. And the II- VI semiconductors are less attractive because the superex- 
change interaction of the doped magnetic ions favors the AFM or spin glass configuration. 
However, in some DMSs based upon the transition metal oxide, the ferromagnetism have 
been observed at or even higher than room temperature, such as 280 K for ZnO^ and 400 K 
for Ti02.^ On the other hand, in the group-IV semiconductor. Park et al7 reported that Mn 
doped Ge (Mnj.Gei_a;) has its Curie temperature up to 116 K and then Cho et al. improved 
it to 285 K.^ So, if it were possible to make room temperature FM Mn^-Sii-a;, the "spin- 
tronic semiconductor" industry would grow up rapidly based upon now mature Si-based 
semiconducting technology and associated facilities. 

The room temperature FM DMSs discovered continuously in experiments brings challenge 
to the theoretical work because their origin of ferromagnetism is still an open question. 
Though some different mechanisms, such as the Ruderman-Kittel-Kasuya-Yoshida (RKKY) 
interaction,^^ double exchange,— double resonance,— the Zener tunneling model,^™ and the 
mean-field theory/^ have been proposed, none could give a conclusive interpretation. The 
strong p — d exchange interaction intermediated by mobile holes is thought as the origin 
of ferromagnetism in III-V compound based DMSs,'- and carrier-induced ferromagnetism 
with the exchange interaction mediated by electrons was considered to be suitable to the 
Co-doped anatase Ti02 system.— While the first-principles calculation on the Mn-doped Ge 
made by Park et al. indicates that the FM order arises from a long-range FM interaction 
competing with a short-range AFM one. More detailed study^^ by Zhao et al. shows that 
Mn2;Gei_x is a RKKY-like FM semiconductor. In this paper, the similar study is extended 
to the Cr (Fe)-doped Ge and Si. Since a Cr (Fe) atom has one electron less (more) than 
a Mn atom and Si is in the same column as Ge in the Periodic Table, comparison between 
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these DMSs allows to investigate variations of magnetic properties with change of dopants 
and host semiconductors. In fact, there have been experiments to be done on the Cr- and 
Fe-doped G o^^i^^ and also the Mn-doped Si.— In the follwoing Section II, we will introduce 
our calculation methods, and our numerical results are shown in the Section III, from which 
some discussions and conclusions are made. 

II. CALCULATION METHOD 

The software package VASP (Ref. 18) has been used in our calculations, which is based 
on a total energy pseudo-potential plane-wave method within the local spin density ap- 
proximation (LSDA). In the calculation, the interaction between ions and electrons is de- 
scribed by the projector-augmented wave method in the generalized gradient approximations 
(GGA).J^ The initial crystal structures of R^Mi-^; (R=Cr, Mn and Fe; M=Ge, Si) are taken 
as 2a X 2a X 2a supercell for x=0. 03125 and 2a x 2a x la for a;=0.0625 with two M atoms 
replaced by two R atoms, among which, the first is at the origin, and the other is put on a 
lattice position farther away from the origin indicated by three digits following N to repre- 
sent its {x,y, z) coordinate in units of a/4.^'* And the lattice constant a is taken as that of 
pure Ge and Si, i.e., 5.658 A and 5.431 A, respectively. In the spin optimization, the initial 
spin configuration in the AFM state is taken as 5 net spin on one R atom, and -5 on the 
other one, and in the FM state, 5 net spins are chosen for both R atoms. The same ground 
state is reached while increasing the value to 8 and 10. An energy of 350 eV is used for 
the plane wave cutoff, and when the energy cutoff is increased to 550 eV, the total energy 
difference between the AFM state and the corresponding FM state changes no more than 
0.05 meV/R. For the Brillouin Zone sampling, we take the same k mesh as that in Ref. 14 
for all the cases. 

III. RESULTS AND DISCUSSION 

The calculated FM total energies and the energy difference between AFM and FM states 
of the Cr-doped Ge and Si are presented in Tables. I and II for x=0. 03125 and 0.0625, 
respectively. Also shown are the averaged magnetic moments on the Cr atom. Obviously, 
for the Cr-doped Ge at both concentrations, the AFM order is energetically more favored 
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than the FM order, which is consistent with the experimental measurement and theoretical 
calculation in Ref. 15. The more holes caused by Cr than Mn in Ge host matrix unexpectedly 
do not enhance the ferromagnetism in this system. And in the Si host matrix, independent of 
the doping concentration, only for the N220 configuration the FM order is a little favorable 
while the others tend to taking the AFM order. But in general the energy differences between 
AFM and FM decreases when x changes from 0.0625 to 0.03125, indicating that perhaps at 
lower enough x, the FM state would be lower than AFM in energy, which still needs more 
experiments to be confirmed. 

Tables. Ill and IV give the results of Mn-doped Ge and Si with a;=0. 03125 and 0.0625, 
respectively. It is clear that our results of Mn^Gei_j; are the same as those in Ref. 14, 
indicating that it is a RKKY-like FM semiconductor. However, it is totally different for 
Mnj.Sii_a; although both Ge and Si are in the same column in the Periodic Table. Inde- 
pendent of the doping concentration, the Mnj,.Sii_2: tends to be FM even when the Mn-Mn 
distance is as long as 9.41 A in the A^444 case except that the A^lll shows AFM order. The 
AFM interaction between Mn ions exists only in a short range (the nearest-neighbor, about 
2.35 A), which competes with the long range (all beyond the nearest neighbor) FM inter- 
action. This mechanism was once thought as the origin of FM order in Mna;Gei_^ system,— 
which now seems to be responsible for the FM order in Mn doped Si. So, it is plausible that 
the FM order would be easier in Si host matrix than in Ge for Mn-doping. When x=0. 03125, 
the lowest FM energy is found in A^440 configuration for Ge and in A^220 for Si, and the 
corresponding energy differences between AFM and FM states {Eafm — Efm) are 103.92 
and 74.88 meV/Mn in these two condigurations. Increasing x to 0.0625, the lowest FM en- 
ergy configuration becomes A^220 for both Ge and Si, and the corresponding Eafm — EpM 
increases to 122.2 and 84.61 meV/Mn, respectively. So, in the same doping concentration, 
the Curie temperature of Mn^-Gei-a; will be higher than that of Mn^-Sii-^^, and in both host 
matrixes, the Curie temperature will increase with Mn doping concentration.— li^^i^iSS 

The results of 3.125% and 6.25% Fe-doped the two group-IV semiconductors are listed in 
Tables. V and VI, respectively. For both concentrations of Fe-doped Ge and Si, the lowest 
energy state are all the A^lll configuration with FM order, and all other configurations have 
much higher energies than the Nlll. So, there would be ferromagnetism in Fe-doped group- 
IV semiconductors, and the recent experimental result shows that Fe^^Gei-i;^^ is a n-type 
FM semiconductor with Curie temperature as high as 233 K while the ferromagnetism in 
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Fe-doped Si still needs further experimental observations. It is found from Table. V, for Ge 
host matrix, all other configurations, except the A^lll and N220, converge to non-magnetic 
state even from the initial FM one, and their energies are lower than that of AFM state. 
In contrast, the N220 favours the AFM state rather than the non-magnetic. However, 
for Si host matrix, all configurations, except only A^lll, favour the non-magnetic, being 
independent of the initial magnetic state. Comparing these results with those at 6.25% in 
Table. VI, we found that the higher Fe concentration will enhance the ferromagnetism. For 
examples, the FM state in all Feo.o625 Geo. 9375 configurations are favorable in energy than the 
AFM. So in the Fe-doping cases, the FM interaction is a very short-range one, and would 
become longer-range by increasing the doping concentration. 

It is generally thought that the d orbital on the R atom is much more hybridized with p 
orbital of Si than with that of Ge, which would cause a smaller local magnetic moment on R 
in Si matrix than that in Ge one.^^ From our calculated magnetic moments on R atoms listed 
in the above mentioned tables it can be seen that at the same concentrations, in general, 
the magnetic moment on R in Ge host matrix is a little larger than that in Si one. We have 
also checked that if the lattice constant of Ge is compressed to that of Si, i. e., increasing 
the mixing of the p orbital of Ge with the d orbital of R atom, the magnetic moment on R 
atom would decreas in all cases. Especially in the case of 6.25% Fe-doping, the decrease of 
Ge lattice constant will cause the Fe atom to be non-magnetic as like that in Feo.0625Sio.9375 
listed in the right part of Table. VI. And vice versa, enlarging Si lattice constant to that of 
Ge causes the A^lll and A^220 of the Fe-doped Si to be in FM order although A^400 and 
N440 reamin to be in the non-maganetic. So, compressing (enlarging) lattice would decrease 
(enhance) the FM order of the TM-doped group-IV semiconductor, which is consistent with 
the theoretical calculation on the Mn-doped diamond, another group-IV element with the 
same crystal structure, in which the absence of ferromagnetism had been predicted.— 

To further investigate the transition metals (TMs) doped group-IV semiconductors, we 
have also studied the electronic structures of these systems, and focused our attention to the 
lowest energy FM states in all the configurations. The obtained total DOS and 3(i-partial 
DOS on R atoms in our concerned systems are shown in Figs. 1 and 2 for a;=0. 03125 and 
0.0625, respectively. Since Cr (Fe) has one electron less (more) than Mn atom, a careful 
analysis of the total DOS of Ct^Mi^x (Fe^^Mi-^;) in both figures shows that the Fermi level is 
a little down-(up-) shifted compared with that of Mn^-Mi.j.. Moreover, from the 3(i-partial 
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DOS and the magnetic moment of R atoms listed in Table. I to VI, we can deduce that one 
absent (redundant) electron in Cr (Fe) compared to Mn is in spin up (down) state, which 
can also be seen from the emergency (narrower) of the gap in spin up (down) channel in 
the case of Cr (Fe) doping compared to that of the Mn doping. These changes of total 
DOS around the Ep in the spin channels will lead to different spin polarized conductivity in 
these systems. As shown obviously in Figs. 1 and 2, the Cr- (Fe-) doped Ge and Si system 
is nearly semiconducting (metal) in the total DOS structure, which is consistent with its 
tendency to be in AFM (FM) groud state because the decreasing (increasing) of the carriers 
mobility would suppress (enhance) the FM order in the system. The 3d partial DOS of the 
R atoms are mostly distributed around in the range from -3 to +2 eV. And in the spin 
up channel, the main peaks of Cr in its valence band are centered at about -1 eV while 
those of Mn and Fe are at -2.3 eV,— which are mostly contributed by and tig. But for 
the spin down part, the occupied peaks of Cr and Fe are mainly located at -0.5 eV and that 
of Mn is at about -1 eV, all of which are mostly composed of t^g, whereas the un-occupied 
state around 0.5 eV for Mn and Fe and that at 1.2 eV for Cr have character. So it is 
known from the above results that in general the crystal field splitting Ag^.^jg is about 0.5, 
1.3 and 1.8 eV and the exchange splitting of Cg orbital is about 2.2, 2.8 and 2.8 eV for Cr, 
Mn and Fe, respectively. Clearly the exchange splitting energy is larger than the crystal 
field splitting energy in these systems, leading to their quite strong spin ordering and little 
changes in host's geometrical structure when doped with TMs.^'^'' And above all, we can see 
that for Mn in Ge (Si), the local magnetic moment is contributed by three electrons in tig 
orbital and the orbital is occupied by one itinerant electron. Since a small crystal field 
splitting, both of them are hybridized with the p orbital of Ge (Si) strongly, which enhances 
the carriers mobility and reduces the magnetic moment on Mn atom. For the Cr, lack of the 
one electron makes this system semiconducting and prefer more to stay in the the AFM 
order. While in the Fe-doped cases, the one more electon occupies the tig orbital, near the 
Fermi level, hybridizing with the p eletrons, which makes the gap in the spin-donw channel 
narrower or disapper, causing the magnetic moment on Fe is smaller than that on Mn. 

According to Anderson's s — d mixing model,— the local magnetic moment on the tran- 
sition metal atom doped into the host matrix is determined by three factors: the on-site 
Coulomb interaction U of d electrons, the mixing of d orbital with the delocalized orbital 
s of the host and the energy difference between the d orbital and the Fermi level of the 
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system. So, we have also done the LSDA+f/ calculations within VASP implement, which 
show that the on-site U will enhance the local magnetic moment on R, especially for the 
6.25% Fe-doped Si. A U=3.0 eV would cause a magnetic moment as large as 3.02 /i^ on Fe 
atom in A^400 configuration and the corresponding FM state is energically lower than the 
AFM one by about 31.67 meV/Fe comparing with the non-magnetic state without U. And 
the total density of states (DOS) of its FM state is nearly half-metal, which is consistent 
with the nearly integer magnetic moment on Fe. The same calculation of Mno.0625Geo.9375 in 
A^220 configuration shows that the system remains to be a half-metal?^ though the magnetic 
moment on Mn is increased to 3.77 fiB- The reason why half-metal feature is kept even in- 
cluding U is due to the strong s — d mixing in these TM-doped group-IV semdiconductors. 
Further, the LSDA+f/ calculation of 6.25% Cr-doped Si in A^220 and A^400 configuration 
shows that the AFM state is still more stable than the FM one in energy. 

IV. CONCLUSION 

In summary, we have studied the transition metal doped group-IV semiconductors 
Ra;Mi_^ (R=Cr, Mn and Fe; M=Si, Ge) by the first principle calculations in LSDA and 
GGA formalism. The obtained results show that there exist different ground states for dif- 
ferent R elements. The higher hole concentration induced by Cr than Mn seems to disfavour 
the FM ordering. On the contrary, Fe-doping makes the group-IV semiconductors more FM 
than Mn-doping. It is also found that enlarging the lattice constant would decrease the 
p — d mixing and be benefit to the FM order, which perhaps is another way to get higher 
Curie temperature magnetic group-IV semiconductor in experiments. These systems seem 
to have different mechanism of FM order although their R elements are in the same row and 
both Si and Ge are in the same column in the Periodic Table, which still needs much more 
experimental and theoretical efforts to make clear. 
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TABLE I: Total energy of the FM phase (Efm) in different configurations counted relatively to 
that with lowest FM energy, and the energy difference between AFM and FM state {Eafm — Efm) 
calculated for Cr-doped Ge and Si at a;=0. 03125. Also shown are the average magnetic moments 
on each Cr atom calculated from initial FM {Mfm) and AFM [Mafm) configuration, respectively. 
Energies are all in unit of meV/Cr and the magnetic moment are in the unit of /Lte/Cr. 





Cro.03125 Geo. 96875 


Cro.o3i25Sio.96875 


System 


Efm 


Eafm-Efm 


Mfm 


Mafm 


Efm 


Eafm-Efm 


Mfm 


Mafm 


iVlll 


397.61 


-387.295 


2.10 


±2.78 


237.265 


-347.695 


1.82 


±2.32 


N22Q 


34.045 


-20.65 


2.54 


±2.68 


14.585 


11.42 


2.22 


±2.27 


iV400 





-18.5 


2.54 


±2.57 





-10.8 


2.22 


±2.22 


N224 


21.255 


-28.375 


2.53 


±2.61 


22.795 


-13.425 


2.20 


±2.23 


iV440 


39.365 


-61.555 


2.48 


±2.72 


17.555 


-23.18 


2.18 


±2.29 


iV444 


28.51 


-31.275 


2.51 


±2.61 


15.61 


-10.98 


2.20 


±2.24 


TABLE II: The same as that in Table. I except at x=0.0625 




Cro.0625 Geo. 9375 


Cro.0625Sio.9375 


System 


Efm 


Eafm-Efm 


Mfm 


Mafm 


Efm 


Eafm-Efm 


Mfm 


Mafm 


mil 


370.58 


-381.27 


2.09 


±2.78 


237.185 


-351.02 


1.84 


±2.33 


N22Q 


10.55 


-36.54 


2.50 


±2.69 


8.72 


5.825 


2.21 


±2.27 


iV400 





-69.78 


2.45 


±2.69 





-32.165 


2.17 


±2.29 


iV440 


12.09 


-07.72 


2.14 


±2.72 


C.005 


-25.025 


2.17 


±2.29 
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TABLE III: The same as that in Table. I, except for the Mn-doped Ge and Si. 







Mno.o3i25Ge 


0.96875 




Mno.o3125Sio.96875 


System 


EpM 


Eafm-Efm 


MpM 


Mafm 


EpM 


Eafm-Efm 


Mfm 


Mafm 


iVlll 


263.27 


-257.395 


3.06 


±3.20 


26.18 


-181.135 


1.18 


±2.69 


iV220 


10.195 


76.425 


3.17 


±3.22 





74.88 


2.83 


±2.48 


iV400 


60.6 


-17.525 


3.16 


±3.21 


136.995 


5.38 


2.81 


±2.81 


N224 


70.355 


-3.91 


3.18 


±3.23 


142.61 


23.735 


2.82 


±2.80 


NUO 





103.92 


3.16 


±3.23 


22.34 


125.145 


2.83 


±2.59 


NAM 


94.295 


-27 


3.16 


±3.25 


153.065 


15.645 


2.81 


±2.83 



TABLE IV: The same as that in Tabic. Ill except at x=0.0625. 





Mno.0625Geo.9375 


Mno.0625Sio.9375 


System 


EpM 


Eafm-Efm 


Mfm 


Mafm 


Efm 


Eafm-Efm 


Mfm 


Mafm 


A^lll 


278.975 


-219.87 


3.03 


±3.18 


158.47 


-258.525 


2.40 


±2.69 


iV220 





122.2 


3.15 


±3.15 





84.61 


2.80 


±2.55 


iV400 


4.93 


104.88 


3.14 


±3.19 


23.275 


109.76 


2.76 


±2.63 




32.655 


98.185 


3.13 


±3.19 


42.54 


108.355 


2.81 


±2.65 



TABLE V: The same as that in Table. I except for Fe-doped Ge and Si. 





Feo.03125 Geo. 96875 


Feo.03125 Sio. 96875 


System 


Efm 


Eafm-Efm 


Mfm 


Mafm 


Efm 


Eafm-Efm 


Mfm 


Mafm 


mil 





170.545 


2.22 


±2.17 





95.33 


1.96 


±1.43 


N22Q 


246.62 


-81.595 


0.001 


±2.47 


129.325 


0.03 


-0.001 


±0.111 


NAm 


184.605 


43.49 


-0.005 


±2.13 


81.725 


-0.01 


0.000 


±0.000 


N22A 


204.575 


39.6 


-0.011 


±2.11 


104.61 





0.000 


±0.001 


iV440 


198.14 


7.895 


-0.006 


±2.62 


86.355 





0.000 


±0.006 


iV444 


211.885 


45.69 


-0.004 


±2.15 


97.65 





0.000 


±0.000 



11 



TABLE VI: The same as that in Table. V except at a;=0.0625. 





Feo.0625Geo.9375 


Feo.0625Sio.9375 


System 


EpM 


Eafm-Efm 


Mfm 


Mafm 


EpM 


Eafm-Efm 


Mfm 


Mafm 


A^lll 





156.04 


2.23 


±2.22 





108.97 


1.98 


±1.51 


iV220 


114.445 


38.485 


2.45 


±2.39 


249.405 


-95.02 


1.57 


±0.21 


Ar400 


97.825 


94.645 


2.50 


±2.06 


101.655 


0.015 


0.004 


±0.02 


Ar440 


129.52 


70.22 


2.55 


±2.05 


106.87 


0.045 


0.02 


±0.02 
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Figure Captions 



Fig. 1. The total DOS (thin hne) and projected 3d partial DOS (bold line) for the lowest 
FM configurations at 3.125% doping concentration. The positive DOS value means that of 
spin up while negative is for spin down. The vertical dashed line represents the Fermi level. 

Fig. 2. The same as Fig. 1 but at 6.25% doping concentration. 
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